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Measurements  have  been  made  of  the  atmospheric  velocity  turbu- 
lence  spectrum  using  a mobile,  CO^,  continuous  wave,  Doppler 
LIDAR  system.  The  LIDAR  was  used  in  a heterodyne  configuration 
with  a 30  cm  output  optics  to  produce  a conical  scan.  Velocity 
data  from  this  scan  were  Fourier  analyzed  to  obtain  the  velocity 
power  density  spectrum.  The  -5/3  power  law  of  Kolmogorov  was 
observed  along  with  variations  in  total  energy  which  could  be 
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correlated  with  storm  activity.  In  addition,  departures  from 
the  Kolmogorov  spectrum  were  observed  under  some  conditions. 
It  is  hoped  that  the  total  energy  and  departure  from  the  -5/3 
power  law  may  prove  useful  in  short-term  meteorological  fore- 
casting. 
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i.  jMKaamatt 

Throughout  modern  time,  man  has  strived  for  accurate 
weather  prediction,  with  a need  for  ever  increasing  amounts  of 
wind  velocity  data.  In  the  early  days  of  meteorology,  wind 
velocities  were  measured  at  widely  separated  points  to  deter- 
mine the  large-scale  features  of  the  wind  velocity  field.  This 
provided  only  a coarse  picture  of  a phenomenon  known  to  have 
small  scale  features  of  interest.  With  the  advent  of  the 
electronic  computer,  it  became  feasible  to  process  the  large 
amounts  of  data  that  newly  developed  laser  systems  could  now 
deliver.  Thus  the  analysis  of  minute  details  of  wind  be- 
havior became  practical. 

A coherent  detection,  CC>2 , continuous  wave  Laser  Doppler 
Velocimeter  (LDV)  is  a useful  instrument  for  these  remote  wind 
velocity  measurements.  The  present  paper  describes  the  use  of 
a LDV  to  determine  the  power  spectral  density  of  a turbulent 
wind  velocity  field,  and  thus  illustrates  the  power  of  the  LDV 
for  collecting,  analyzing,  and  displaying  large  amounts  of  wind 
velocity  data. 

The  success  of  the  system  is  due  to  a combination  of  the 
small  scale  remote  measurement  capability  of  the  laser  with  the 
numerical  analysis  capability  of  a computer,  which  can  reduce 
large  quantities  of  data  to  meaningful  displays.  A sample  of 
several  hundred  thousand  wind  velocity  measurements  is  reduced 
to  an  energy  spectral  density  curve  with  easily  recognizable 
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features.  This  system  has  been  used  to  verify  the  -5/3  power 
law  of  Kolmogorov  at  spatial  frequencies  from  1/30  to  1/2  km-1. 

This  system  is  the  result  of  more  than  ten  years  of  re- 
search and  development  of  CO,  laser  systems.  The  first 
measurement  of  the  Doppler  shift  in  CO^  radiation  backseattered 
from  the  atmosphere  was  made  in  1964.  By  1974,  sophisticated 
systems,  using  computers  for  real  time  processing,  had  been 
used  for  measuring  the  behavior  of  aircraft  wakes  and  dust 
devils.  Atmospheric  measurements  have  been  made  at  ranges  up 
to  1 km  with  CW,  20  watt  systems  and  8 kn  with  25  mJ  pulsed 
systems . 

The  coherent  COn  LDV  is  ideally  suited  for  measurement  of 
the  line-of-sight  velocity  component  as  a function  of  range 
and  time.  Measurement  of  a three  dimensional  velocity  vector 
requires  at  least  three  independent  measurements.  These  may  be 
made  directly  by  focussing  three  systems  into  the  same  volume 
of  space  from  different  angles,  or  may  be  inferred  either  from 
successive  measurements  at  a given  point  in  space  with  a moving 
system,  or  from  closely  timed  measurements  at  slightly  differ- 
ent points  in  space  using  a fixed  system.  The  moving  system 
approach  assumes  that  the  velocity  vector  does  not  change  in 
the  measurement  time,  while  the  fixed  system  approach  assumes 
that  the  velocity  vector  does  not  vary  over  the  space  required 
for  the  measurement.  The  present  system  uses  the  fixed  system 
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approach  by  generating  a conical  scan  about  a vertical  axis. 

A 30°  cone  half  angle  is  used  so  the  diameter  of  the  circle 
described  in  space  is  equal  to  the  slant  range  to  the  focus 
of  the  transmitted  beam.  The  velocity  vector  at  the  center  of 
the  circle  is  inferred  from  the  measurements  of  the  line-of- 
sight  velocity  around  its  circumference. 

This  paper  begins  with  a brief  description  of  the  assumed 
wind  velocity  field  and  the  required  measurements  and  calcula- 
tions. The  principles  of  obtaining  velocity  measurements  with 
an  LDV  using  a conical  scan,  and  the  principles  of  using  this 
velocity  to  obtain  the  turbulent  velocity  energy  spectral  den- 
sity are  described.  Then  the  system  is  described  in  detail, 
followed  by  the  processing  methods,  results,  and  conclusions. 

It  will  be  shown  that  the  LDV  provides  a useful  tool  for  mea- 
surement of  mesoscale  atmospheric  behavior  with  a capability 
for  obtaining,  processing,  and  analyzing  large  quantities  of 
data  concerning  wind  velocity  fields. 

2.  THE  TURBULENT  KINETIC  ENERGY  DENSITY  SPECTRUM 

The  wind  velocity  is  a complicated  function  which  varies 
randomly  with  position  and  time.  Thus,  even  if  the  wind  veloc- 
ity is  accurately  measured  at  a large  number  of  points,  the  re- 
sults may  not  be  easily  interpreted.  Some  method  of  simpli- 
fying the  data  to  present  meaningful  results  in  usable  form  is 
required.  One  approach  to  this  problem  is  Fourier  analysis. 

The  behavior  of  the  wind  velocity  at  a given  time  may  be 
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expressed  by  the  vector  velocity  9(r)  as  a function  of  position, 

at  . 

or,  equivalently  by  its  spectrum  V (k)  where  k is  defined  as  a 
spatial  frequency.  These  variables  are  a Fourier  transform 
pair: 


V(k) 


(2h) 


3/2 


///*<;>  d3r 


While  these  two  representations  contain  identical  informa- 
tion, it  is  reasonable  to  expect  that  one  or  the  other  may 
present  the  information  in  a way  which  is  more  understandable 
and  more  directly  related  to  the  physics  of  the  problem.  For 
the  present  purpose,  it  is  the  velocity  spectrum  which  proves 
to  be  most  useful.  Representation  of  the  wind  velocity  in 
terms  of  its  spectrum  is  equivalent  to  stating  that  the  veloc- 
ity field  consists  of  many  sinusoidal  waves  having  different 
amplitudes  and  phases.  The  square  of  the  amplitude  | V (k) | “ has 
units  of  jenergy  mass  * (spatial  frequency)  "|J,  and  represents 
the  amount  of  Kinetic  energy  per  unit  mass  of  atmosphere  per 
unit  volume  of  frequency  space. 

It  should  be  evident  that  this  would  be  useful  in  studying 

the  distribution  of  Kinetic  energy  in  the  atmosphere.  Theoret- 

12  3 

ical  treatment  of  the  problem  by  Kolmogorov  ' ' has  indicated 
a -5/3  power  law  should  exist  for  this  function; 


I 


4 


SS& 


RAYTHEON 

CQUIPMKNT 


COMPANY 

o i v i • i o n 


E • | v(k)  I 2 


A£2/3k‘5/3 


where  A is  a constant  of  magnitude  near  h and  e is  the  rate  of 
dissipation  of  turbulent  energy  per  unit  mass.  Kinetic  energy 
normally  enters  the  atmosphere  at  low  spatial  frequencies  (large 
wavelengths) , as  a result  of  interaction  of  the  atmosphere  with 
large  features  of  the  landscape,  pressure  gradients  and  other 
large-scale  causes.  Kinetic  energy  is  removed  by  conversion  to 
thermal  energy  most  easily  at  small  wavelengths.  While  these 
processes  occur,  the  atmosphere  redistributes  energy  throughout 
the  spectrum  in  an  attempt  to  maintain  the  -5/3  power  law.  When 
a change  in  the  energy  balance  occurs,  the  -5/3  power  law  should 
be  temporarily  upset  at  certain  wavelengths  until  the  equilib- 
rium is  reestablished.  Thus,  one  could  speculate  that  measure- 
ment of  the  energy  spectral  density  should  yield  information  on 
impending  weather  changes,  by  indicating  change  in  the  energy 
balance  of  the  atmosphere.  In  particular,  an  increase  in  the 
energy  at  long  wavelengths,  above  that  anticipated  by  the  -5/3 
law,  could  indicate  that  the  total  energy  in  the  atmosphere  is 
increasing.  Such  an  increase  would  most  likely  be  followed  by 
a storm  in  which  the  energy  could  redistribute  itself  according 
to  the  -5/3  power  law.  Thus,  the  energy  spectrum  of  turbulent 
wind  velocity  is  an  important  variable  in  meteorological 
research. 


It  is  often  convenient  to  use  the  autocorrelation  function 
of  a random  variable  as  a way  of  characterizing  the  size  scale 
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over  which  the  variable  changes.  This  function  is  the  Fourier 
transform  of  the  energy  density  spectrum  and  may  readily  be 
obtained  from  experimental  data.  The  autocorrelation  function 
is  given  by 

r(Ar)  = <V(r)  V(r-tr)> 

where  the  brackets  denote  an  ensemble  average.  If  the  atmos- 
phere is  homogeneous,  the  ensemble  average  may  be  replaced  by  a 
spatial  average  over  a sample  of  space.  This  provides  a simple 
way  of  obtaining  the  correlation  function  and  thus  the  energy 
spectrum.  The  correlation  function  in  one  dimension  is  ob- 
tained by  averaging  over  x as 

I*  (Ax)  = (X  -x~T  f*2  ^(x)*^(x  ‘ Ax)  631 

xi 

The  energy  spectrum  is  then  simply  the  Fourier  transform  of  the 
function. 

Finally,  it  has  been  postulated  that  it  is  possible  to 

measure  r(Ax)  at  a single  point  in  space  by  making  a measure- 

4 

ment  over  a period  of  time  and  invoking  Taylor's  Hypothesis. 
This  hypothesis  states  that  the  temporal  and  spatial  correla- 
tion functions  are  equivalent  if  the  time  is  multiplied  by  the 
mean  wind  velocity; 

<V (x , t)  V (x , t+T)>  = <V (x , t)  V (x-ut  , t) > 
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This  is  often  referred  to  as  the  assumption  of  frozen  turbu- 
lence. That  is,  the  turbulent  wind  field  as  viewed  by  an  ob- 
server moving  at  the  mean  wind  velocity  changes  more  slowly 
than  that  measured  by  a stationary  observer.  The  gusts,  viewed 
from  a stationary  position,  are  thus  slowly  changing  eddies 
being  transported  by  the  mean  wind.  This  assumption  allows  a 
particularly  simple  measurement  of  j V | as  a function  of  time, 
at  one  location,  to  be  used  to  construct  the  correlation  func- 
tion, r(x) . Furthermore,  for  isotropic  turbulence,  this  leads 
directly  to  the  three  dimensional  correlation  function  r(r),  and 
thus  to  the  energy  density  spectrum  E (k) . 

3.  PRINCIPLES  OF  OPERATION  OF  THE  SYSTEM 


The  principles  of  operation  of  the  Coherent  CO^  LDV  in- 
clude three  areas  of  interest:  the  first  is  the  use  of  the 
Doppler  shift  to  obtain  a line-of-sight  velocity  measurement, 
the  second  is  interpretation  of  a series  of  different  line-of- 
sight  velocity  measurements  to  obtain  a three  dimensional  veloc- 
ity vector,  and  the  third  is  the  use  of  a series  of  vector 
velocities  obtained  at  different  times  to  infer  the  turbulent 
velocity  power  density  spectrum.  In  this  section,  each  of  these 
will  be  discussed  in  turn. 


The  measurement  of  a single  component  of  velocity  by  the 
coherent  CC>2  LDV  may  be  understood  with  the  aid  of  Figure  1. 

A laser  beam  is  transmitted  through  an  interferometer,  which 
passes  most  of  the  energy  to  the  target,  and  diverts  a small 

amount  to  the  detector  as  shown  in  Figure  la.  The  transmitted 
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A)  PATHS  OF  TRANSMITTED  BEAM 


B)  PATHS  OF  TRANSMITTED 

AND  RECEIVED  BEAMS  EOA-106 

Pig.  1.  Doppler  LIDAR  Transmission 
and  Reception. 


WIND  VELOCITY  FIELD 


EOA-107 


Pig.  2.  Velocity  Component 

Sensed  by  the  LIDAR. 


REFERENCE 

BEAM 


/~\ 


SIGNAL 

BEAM 


V<I  » — 


EOA-108 


Fig.  3.  Spectral  Distribution  of 
Paver  on  the  Detector. 


FREQUENCY(MHi) 


VELOCITY  (METERS  PER  SECONO) 

EOA-109 

Fig.  4.  Typical  Lidar  Doppler 
Signal  Spectrum. 
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beam  is  usually  passed  through  a telescope  to  the  target,  which, 
in  the  present  case,  consists  of  naturally  suspended  aerosols 
in  the  atmosphere.  Some  of  the  transmitted  energy  is  back- 
scattered  by  the  target  particles,  and  enters  the  detector  by 
the  path  shown  in  Figure  lb.  According  to  the  Doppler  princi- 
ple this  returning  energy  has  been  shifted  in  frequency  as 
illustrated  in  Figure  3.  When  these  beams  are  superimposed, 
the  result  contains  a contribution  at  the  difference  frequency, 


The  difference  frequency  is  the  Doppler  frequency,  which  is 
proportional  to  the  magnitude  of  the  target  velocity  com- 
ponent, V|j . A sample  spectrum  in  Figure  4 shows  the  intensity 
as  a function  of  frequency  for  an  actual  LDV  atmospheric  return. 
It  may  be  seen  that  there  is  a peak  at  .7  MHz  corresponding  to  a 
velocity  of  about  3.3  meters  per  second.  It  will  be  noticed  that 
the  sign  of  the  velocity  has  been  lost  at  this  point. 


The  upward  looking  VAD  (velocity-azimuth  display)  scan  is 
useful  for  measuring  three  components  of  the  wind  velocity  at 
a point  above  the  system.  The  scan  configuration  is  as  shown  in 
Figure  5.  The  focussed  laser  beam  is  made  to  scan  a cone  about 
a vertical  axis  with  the  focus  at  a range,  R,  giving  the  desired 
altitude,  h.  In  a uniform  wind  field,  the  velocity  along  the 
line-of-sight  for  a vertical-axis  conical  scan  is: 
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Vl|  * u cos  0 sin  d + v sin  0 sin  d + w cos  d 
where  u,  v,  and  w are  the  x,  y,  and  z components  of  velocity, 
0 is  the  azimuth  angle  from  the  x axis,  and 
d is  the  cone  half  angle  from  the  z axis. 

It  is  easy  to  show  that  a least  squares  fit  to  a set  of  data 
consisting  of  velocities  and  azimuths,  and  0^,  such  that 


is  a minimum,  produces  a set  of  equations  expressible  in  matrix 
form  as: 

S = MY 


2 

cos  0 


sin  0 cos  0 
cos  0 


sin  0 cos  0 
2 

sin  0 
sin  0 


cos  0 

sin  0 
1 


u sin  d 

v sin  d 
w cos  d 


) 


where  the  bar  over  a quantity  denotes  the  average  as  in  this 

example:  N 

V cos  0 = •“  ^ cos  0i 

i=l 

The  problem  of  calculating  u,  v,  and  w from  the  and  0^  is 
thus  reduced  to  the  simple  task  of  inverting  the  matrix  JJ: 


2 


§. 
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If  the  points  are  uniformly  spaced  over  rr  or  2n  radians,  the 
matrix  becomes 

(500 
0 .5  0 
0 0 1 

and  the  data  processing  involves  only  calculation  of  the  matrix 
V and  a few  simple  mathematical  operations. 

A sample  of  the  data  for  a typical  VAD  scan  is  shown  in 
Figure  6.  The  quantized  curve  is  the  raw  data  with  its  sign 
reconstituted  by  estimating  the  wind  direction  using  the  high- 
est velocity  and  associated  angle  with  the  assumption  of  zero 
vertical  wind.  The  smooth  curve  shows  a least  square  fit  to  a 
sine  wave  which  was  performed  in  real-time,  using  a simple 
algorithm.  Figure  7 shows  a set  of  u,  v,  and  w components  as 
functions  of  time  from  data  collected  with  a VAD  scan.  It  is 
obvious  in  Figure  6 that  some  errors  occurred  in  the  sign  re- 
constitution process.  These  may  be  eliminated  by  an  iterative 
process  which  uses  the  calculated  u,  v,  and  w to  determine 
corrected  zero  crossings  followed  by  a new  least  squares  fit. 
The  process  will  reach  a final  value  usually  within  two  or 
three  iterations  and  will  seldom  change  u,  v,  or  w by  more  than 
2%. 

It  should  be  noted  at  this  point  that  the  data  includes 
360  data  points,  from  which  360  independent  parameters  may  be 
obtained.  This  may  be  considered  in  the  following  way.'’  If 
the  wind  field  in  the  plane  of  the  scan  can  be  expressed  as  a 
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power  series 

<*>  00 

V (x,y>  - V (o.o)  + £ *xnX"X  Vy" 

n=l  m=l 

where  the  origin  of  coordinates  is  at  the  scan  center  and  a 

and  a are  unknown  constants,  then  the  Fourier  series  of  V ( © ) 
yn 

gives  information  about  V (o,o)  from  the  zeroth  and  first 

terms,  and  about  a and  a from  the  nth  term.  Calculation 

xn  yn 

of  V(o,o)  = ux  + dy*  + wz'  is  the  same  as  in  the  above  least 
square  fit  algorithm. 

This  has  two  interesting  results:  first,  it  states  that 
it  is  possible  to  obtain  some  information  on  the  wind  veloc- 
ity field  inside  the  circle,  and  second,  it  states  that  it  is 
not  necessary  to  make  the  assumption  of  a uniform  wind  field 
over  the  scan  in  order  to  obtain  a good  estimate  of  the  veloc- 
ity vector  at  the  scan  center.  Thus,  the  system  actually 
measures  an  average  velocity  over  a region  smaller  than  the 
circle. 


•J  vM  W * OJ  W*JWVrl  ; 


The  coherent  CC>2  LDV  is  a laser  system  capable  of 
measuring  line-of-sight  atmospheric  velocity  using  the  Doppler 
shift  as  described  in  the  previous  section.  The  optical  com- 
ponents transmit  the  laser  radiation,  receive  it,  and  recombine 
it  with  the  reference  beam  to  generate  the  frequency  shift. 

The  LDV  processor  measures  this  frequency  shift,  digitizes  and 
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inputs  the  data,  including  the  spatial  coordinates  of  the 
point  at  which  the  velocity  was  measured,  to  a mini  computer. 
Through  the  use  of  the  mini  computer,  the  system  is  capable  of 
operation  at  a remote  site  with  a guick  look  capablity  for 
system  performance  and  data  quality  evaluation.  The  mini  com- 
puter performs  real  time  processing  to  prepare  the  data  for 
output  to  a large  computer  for  the  final  processing.  A block 
diagram  of  this  system  is  shown  in  Figure  8 where  each  of  the 
blocks  are  described  below: 

The  laser  is  a Raytheon  model  LS-10A  with  a nominal  out- 

2 

put  power  of  17  watts.  The  output  beam  is  TEMqo  with  the  1/e 
intensity  points  at  a diameter  of  about  7 mm.  The  horizontal 
polarization  is  determined  by  Zinc  Selenide  (ZnSe)  Brewster 
windows  on  the  ends  of  the  discharge  tube. 

In  the  interferometer,  the  reference  beam  is  obtained 
from  the  transmitter  and  combined  with  the  reflected  radiation. 
The  telescope  is  an  f/8  Cassegrainian  with  a one  foot  primary. 
It  serves  the  purpose  of  expanding  the  beam  and  controlling  the 
range  to  focus  by  means  of  a moving  secondary.  The  detector 
is  a liquid  nitrogen  cooled  (77K)  Lead  Tin  Telluride  chip  of 
about  250  nin  width.  It  is  mounted  in  a dewar  capable  of 
holding  liquid  nitrogen  for  8 hours.  The  reference  beam  power 
of  a few  milliwatts  is  sufficient  to  provide  shot  noise  limited 
operation  of  the  detector  over  a bandwidth  greater  than  foO  MHz. 
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The  spectrum  analyzer  and  spectrum  processor  consist  of 
a SAW  delay  line  spectrum  analyzer  followed  by  a digital  pro- 
cessor which  provides  for  temporal  and  spectral  integration, 
amplitude,  low  frequency  and  frequency  width  thresholds,  linear 
or  logarithmic  output  ranges  and  digital  output  of  the  selected 
spectrum  and  parameters.  The  SAW  delay  line  is  ideally  suited 
for  spectral  analysis  because  of  its  high  duty  cycle  frequency 
measurements.  The  processor  bandwidth  is  10  MHz  with  a reso- 
lution of  100  kHz.  The  amplitude  has  256  discrete  levels. 

Parameters  which  are  available  to  the  computer  are: 

Ip^,  the  highest  amplitude  in  the  thresholded 
spectrum, 

Vm>,  the  frequency  (velocity)  associated  with  Ip^, 

V the  highest  velocity  in  the  threshold  spectrum, 
and 

N,  the  number  of  velocity  filters  above  threshold. 
These  are  illustrated  in  Figure  9. 

The  mini  computer  is  a Raytheon  706  with  a cycle  time  of 
1 usee  and  16K  of  core  memory.  It  is  equipped  with  paper  tape 
facilities,  a Tektronix  4010-1  display,  » hard  copy  unit,  a 
time  code  generator,  and  an  interface  to  collect  data  from  the 
system. 

For  final-processing,  the  CDC-6700  computer  installation 
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accessed  via  phone  lines  to  perform  more  sophisticated  data 

1 

processing  with  Fortran  programs  and  plotting  routines  for 
display  of  the  processed  data. 

The  scanner  consists  of  a tilted  mirror  rotating  about  a 
vertical  axis.  The  beam  enters  from  above  using  a bending 
mirror,  and  then  is  reflected  from  the  scan  mirror  at  a 30° 
angle  to  the  vertical.  As  the  scanner  rotates,  the  focussed 
spot  describes  a circle  in  a horizontal  plane  above,  and 
centered  on,  the  system. 

Figure  10  shows  a view  of  the  van  in  which  the  system  is 
mounted . 


The  data  obtained  from  the  coherent  C02  LDV  is  processed 
using  a real-time  correlation  routine  on  a Raytheon  706  mini 
computer  followed  by  an  FFT  performed  on  a CDC-6700 . The  FFT 
routine  includes  numerous  displays  for  evaluation  of  the  data. 

An  overall  flow  of  the  processing  is  shown  in  Figure  11  in 
two  parts.  The  top  part  is  the  processing  carried  out  by  the 
Raytheon  706  mini  computer  installed  with  the  system.  The  lower 
part  is  the  processing  accomplished  by  the  CDC-6700  at  the 
Raytheon  installation  in  Bedford,  Mass. 

Use  of  the  two  computers  permits  maximum  flexibility  in 
the  final  processing  algorithms  while  retaining  a quick-look 
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Figure  11.  Flow  Diagram  of  Data  Processing 
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capability  in  the  real  time  data  processing  phase.  The  two 
processing  operations  are  linked  by  one  of  two  data  transfer 
techniques:  by  physically  carrying  the  paper  tape  containing 

the  correlation  functions  to  a terminal  that  could  read  it  into 
the  CDC-6700 , or  by  feeding  the  information  directly  to  the  6700 
through  a telephone  link  connected  to  the  706.  The  following 
paragraphs  discuss  separately  the  processing  algorithms  used 
by  the  two  computers. 

The  velocity  information  V(0)  received  from  a conical  scan 
in  the  presence  of  a uniform  horizontal  wind  field  would  be 
sinusoidal  with  a peak  occurring  when  the  azimuth  of  the  beam  is 
upwind,  a trough  when  it  is  downwind  and  zeros  when  it  is  per- 
pendicular to  the  wind  field.  The  presence  of  a vertical  wind 
causes  a constant  offset  of  this  sine  wave.  The  present  LDV, 
as  discussed  earlier,  cannot  detect  the  sign  of  the  line-of- 
sight  velocity  and  this  ambiguity  results  in  a rectification 
of  the  sinusoid.  A velocity  vs.  azimuth  graph  for  a uniform 
wind  field  with  an  updraft  present  is  shown  in  Figure  12,  as  it 
would  be  seen  by  the  LDV  after  completing  one  360°  scan  around 
the  cone.  The  information  to  be  obtained  is  the  magnitude  of 
the  local  velocity  averaged  over  a scan.  The  technique  for 
obtaining  this  velocity  was  outlined  in  the  principles  of 
operation.  It  is  implemented  in  the  following  way  (see 
Figure  13) . First  360  points  are  collected  and  the  largest 
value  is  found.  It  is  then  assumed  that  any  points  outside  of 
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Figure  12.  Idealized  Raw  Data 
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Figure  13.  Velocity  Calculation 
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180°  bin  centered  on  the  largest  value  were  originally  of  oppo- 
site sign.  In  the  reconstruction  step  all  of  these  rectified 
points  are  made  negative.  The  data  is  now  in  the  form  of  a 
rough  sinusoid  which  is  least  square  fitted  to  the  function: 


Vj|  (6^)  = u cos  0 sin  6 + v sin  0 sin  6 + w cos  6 
Next,  the  magnitude  of  the  velocity  is  calculated; 

| V)  = / u2  + v2  + w2  * 

Then,  the  power  spectral  density  is  calculated.  This 
could  be  accomplished  by  squaring  the  magnitude  of  the  Fourier 
transform  of  the  velocity  data,  but  that  would  require  that 
every  data  point  be  saved  so  that  the  transform  could  be  taken 
at  the  end  of  the  run.  During  a typical  run,  this  would  amount 
to  many  thousands  of  points.  A method  employing  more  economi- 
cal usage  of  computer  memory  is  to  take  Fourier  transform  of  the 
autocorrelation  of  the  velocity  data  as  described  earlier.  A 
1024  point  autocorrelation  is  employed. 


The  correlation  function 


is  computed  by  the  equation: 


r(J) 


N, 


max 


N 


max 


L iv  i 

N*»l 


where  J is  an  index  from  1 to  1024 

N is  the  total  number  of  data  points 

luclX 

The  process,  shown  in  Figure  14,  is  as  follows:  a 1024  word 
array  is  filled  up  with  the  first  1024  velocities,  |V^|  through 


j 

I 


I 

( 


I 


24 


I 


O I V I • t o N I 


COttFlATION 


Figure  14.  Correlation  Process 
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| vio24 I ' and  ^or  eac^  additional  data  point  the  appropriate 
products  are  added  to  the  sums  in  the  above  equation.  The  run 
is  terminated  by  the  system  operator  toggling  a switch  on  the 
706  computer  when  enough  data  has  been  taken.  The  correlation 
function  is  then  output  from  the  computer  to  the  hard  copy  unit 
and  the  large  computer. 

The  final  data  processing  is  performed  on  Raytheon's 
CDC-6700  computer  to  determine  the  energy  density  spectrum. 

The  first  function  of  the  processing  program  is  to  check  the 
data  format  to  determine  if  any  errors  occurred  during  the  data 
transmission  or  paper  tape  reading  process.  An  opportunity  is 
provided  for  correcting  errors  at  this  point.  Next,  the 
program  subtracts  the  average  value  of  the  correlation  function 
to  remove  the  large  DC  contribution  to  the  transform.  At  this 
point  the  data  consists  of  1,024  values  of  the  correlation 
function  correctly  sequenced  in  time,  with  the  spacing  between 
data  points  corresponding  to  the  time  required  for  one  scan. 

These  values  are  plotted.  The  next  function  of  the  processing 
program  is  Hamming  weighting.  Whenever  a finite  data  sample  is 
collected,  the  Fourier  transform  is  equal  to  the  Fourier  trans- 
form of  the  desired  function  convolved  with  the  Fourier  transform 
of  a function  which  is  unity  for  the  sample  interval  and  zero 
elsewhere.  The  side  lobes  of  this  function  contribute  to  a 
noise  floor  in  the  spectrum.  This  effect  is  reduced  by  multi- 
plying the  data  by  a function  which  has  a Fourier  transform 
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with  a broader  peak,  but  much  lower  side  lobes.  The  function 
used  in  the  present  case  is  the  following: 

CO.2  (*&*-  ) * -OB 

This  function  results  in  a signal-to-noise  ratio  of  approximate- 
ly 30  dB  in  the  spectrum.  After  the  Hamming  weighting  is  per- 
formed, a plot  of  the  data  is  again  produced  to  show  the  effects 
of  the  weighting.  The  correlation  function  is  known  to  be  sym- 
metric about  T * 0,  and  this  results  in  the  equivalent  of  a 
2048  point  transform.  The  Fourier  transform  is  accomplished 
using  a simple  fast  Fourier  transform  (FFT)  algorithm.  After 
the  Fourier  transform  is  taken,  the  data  consists  of  spectral 
data  from  zero  frequency  to  a frequency  corresponding  to  one- 
half  the  inverse  of  the  sample  time.  Since  the  correlation 
function  is  symmetric,  the  power  spectral  density  is  twice 
the  real  part  of  the  Fourier  component.  These  values  are 
printed  in  tabular  form  and  are  also  displayed  on  a logarithmic 
plot  with  appropriate  scales.  The  frequency  scale  is  calibrated 
in  units  of  (uT)”*  where  T is  the  sample  separation.  Also  shown 
on  this  plot  is  the  resolution  size  based  on  the  size  of  the 
circle  which  is  scanned  by  the  system.  Displaying  the  data  in 
a logarithmic  format  makes  it  possible  to  determine  the  ex- 
ponent of  the  spectral  density,  the  amount  of  energy  present  in 
the  spectrum,  and  any  deviations  from  the  expected  power  law. 
Thus,  the  final  processing  algorithm  provides  the  data  in  a use- 
ful format  for  interpretation  as  it  relates  to  meteorological 
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data  associated  with  the  experiment. 


6.  &5SI2LIS 

Thirty  tests  were  conducted  during  the  time  period  from 
August  1976  through  March  1977.  Of  these,  17  used  a simpler 
velocity  processing  algorithm  where  the  highest  velocity  in  a 
scan  is  used  as  |V| , while  13  used  the  least  square  fit  to  a 
sine  wave.  These  13  were  collected  during  February  and  March 
of  1977.  Seven  runs  in  March  were  correlated  with  daily 
meteorological  data  obtained  at  Boston,  Massachusetts,  about 
32  km  east  of  the  test  site.  For  each  of  these  runs,  the 
wind  speed  was  averaged  and  a standard  deviation  obtained  to 
supplement  the  correlation  function  and  power  spectral  density 
plots.  The  altitude  of  these  measurements  was  500  meters.  A 
sample  of  these  plots  is  shown  in  Figure  15.  The  correlation 

2 

function  has  it  average  value  removed,  and  has  units  of  knots  . 

The  peaks  of  the  power  spectral  density  plot  show  a -5/3  slope 

2 3 

yielding  an  € of  .003  cm  /sec  . The  weather  at  the  time  of  the 
run  was  clear  and  warm  following  a storm  the  previous  day  which 
produced  0.6  inches  of  rain. 

Figure  16  summarizes  the  results  of  the  data  collected 

during  March  with  the  meteorological  data.  It  may  be  seen  that 

-4  2-3 

c varies  from  near  10  to  .04  cm  sec  . High  values  of  e 
tend  to  correlate  with  low  barometer  readings  and  precipitation. 
A phenomenon  which  is  observed  on  several  occasions  is  a spec- 
trum with  two  regions  of  nearly  -5/3  slope  having  different 
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values  of  e . This  is  illustrated  in  Figure  17  which  is  from 
run  27.  This  run  was  conducted  on  29  March?  a cloudy  day  with 
showers.  A similar  phenomenon  was  observed  in  run  36  on  31 
March  under  similar  conditions  as  shown  in  Figure  18. 

It  is  worthwhile  to  evaluate  some  of  the  data  collected 
using  the  simpler  algorithm  since  it  includes  two  runs  during  a 
hurricane.  These  runs  were  collected  on  August  1976  while 
Hurricane  Belle  was  between  Long  Island,  New  York  and  the 
Massachusetts-Vermont  state  line.  Unfortunately,  the  simplici- 
ty of  the  velocity  algorithm  leads  to  some  erroneous  results, 
including  a slope  less  steep  than  -5/3.  However,  the  discontin- 
uities described  above  are  clearly  visible  in  Figure  19  with  the 
most  structure  seen  in  run  1.  After  landfall,  the  hurricane  de- 
creased in  intensity  and  the  structure  in  the  later  run  (2)  was 
greatly  reduced.  For  comparison,  run  3 was  made  several  days 
later  in  clear  weather,  and  shews  very  little  structure. 

In  summary,  several  good  sets  of  data  have  been  obtained. 
The  Kolmogorov,  -5/3  power  law  is  observed  to  be  met  fairly 
well  on  most  runs.  The  total  turbulent  energy  increases 
during  storms,  and  the  spectrum  often  shows  multiple  regions  of 
-5/3  slope  which  appear  to  be  correlated  with  storms.  While 
considerably  more  data  must  be  collected  and  analyzed  to  de- 
termine the  usefulness  of  the  energy  density  spectrum  in  fore- 
casting, the  coherent,  C02  LDV  has  been  demonstrated  to  be  a 
powerful  tool  for  making  meteorological  measurements. 
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WEATHER  / cloudier, 
may  shower,  low  60s 

BOSTON  AREA  - Becoming  pertly  cloudy,  chance  of  a 
lew  showers.  highs  in  the  low  80*  (H  to  U Celsius), 
southwesterly  winds  15  to  35  mph  Tonight,  fate.  Iowa  in  tow 
40s  (4  to  7 C)  Tomorrow.  partly  tunny,  highs  around  70  (21 
C>.  Probability  a t precipitation  40  peroont  today  and  IB 
percent  tonieht 
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